Quantum spectrum as a time series : Fluctuation measures 



O 
O 

(N 



o 

m 

Q 

U 
d 



> 

in 
m 
O 
oo 
O 
in 
o 



> 
X 



M. S. Santhanam, Jayendra N. Bandyopadhyay and Dilip Angom 
Physical Research Laboratory, Navrangpura, Ahmedabad 380 009, India. 

The fluctuations in the quantum spectrum could be treated like a time series. In this framework, 
we explore the statistical self-similarity in the quantum spectrum using the detrended fluctuation 
analysis (DFA) and random matrix theory (RMT). We calculate the Hausdorff measure for the 
spectra of atoms and Gaussian ensembles and study their self-affine properties. We show that DFA 
is equivalent to A3 statistics of RMT, unifying two different approaches. We exploit this connection 
to obtain theoretical estimates for the Hausdorff measure. 
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The fluctuations in physical systems carry important 
information about the system. In the context of quan- 
tum systems, the spectral fluctuations reveal whether the 
corresponding classical dynamics is regular or chaotic or 
a mixture of both [l| . It is well known that for the clas- 
sically regular systems, the quantum eigenvalues are un- 
corrected while for the chaotic systems, the eigenvalues 
tend to display a certain degree of correlation, which is 
dictated by the symmetry properties of the system such 
as the presence or absence of rotational and time reversal 
invariance. Typically, the spectrum of high dimensional 
quantum systems like the atoms, molecules and nuclei 
belong to the latter class. This equivalence between the 
nature of classical dynamics and spectral fluctuations in 
the corresponding quantized system is generally through 
an analogy with random matrix ensembles [lj and has 
been verified in many simulations and experiments 2]. 

Recently a different approach has been suggested. It 
is possible to consider suitably transformed eigenvalues 
of a quantum system as a time series. Then, using time 
series analysis methods, it was shown that the quantum 
systems display /~ 7 noise where 1 < 7 < 2 depends on 
the classical dynamics in the system [3|. In particular, 
for the chaotic systems the level fluctuations display 1// 
noise. So do the levels of complex nuclei In contrast, 
for the regular systems we get l// 2 noise. Hence, the 
spectral fluctuations can be characterized without refer- 
ence to random matrix models. 

This is reminiscent of following types of time series. 
In the Brownian motion or random walk time series 
Xt, t = 1, 2, . . . , N, the successive increments of the se- 
ries are uncorrelated. The variance of this process is 
var(a;) oc t 2H , (H — 1/2), where H is the Hausdorff mea- 
sure. In this case, the power spectrum gives rise to l// 2 
noise j|. Thus, regular classical systems with uncorre- 
lated quantum levels are analogous to a series of Brow- 
nian motion path. On the other hand, one could also 
imagine a time series with 1// noise. There are many 
mechanisms and processes that produce 1// noise In 
a chaotic system, the eigenvalues of the corresponding 
quantum system contains certain degree of correlation 
and it displays 1// noise. In this case, H — and corre- 
sponds to variance being time- independent. In general, 
for < H < 1,7 = 2H + 1. In terms of the exponent H, 
our interest is in the regular (H = 1/2) and chaotic limits 



(H = 0) for the fluctuations of the spectral 'time series'. 
This region < H < 1/2 corresponds to anti-persistent 
time series, i.e, the one which has opposing trends at suc- 
cessive time steps. Visually an anti-persistent time series 
presents a very rough profile as compared to H = 1/2 se- 
ries. In fact, H is also used to characterize the roughness 
of surfaces and profiles. 

Such an anti-persistent series is also a self-affine frac- 
tal and displays statistical self-similarity that scale differ- 
ently in different directions. Mathematically, it implies, 
y(Ht) k t~ H y(t) where the symbol ~ represents statisti- 
cal equality. Note that the power spectrum is a measure 
of the strength of frequencies in Fourier space but the 
fluctuation function and the Hausdorff measure, taken 
together, indicate the structure of the spectrum at vari- 
ous spectral scales. This structure is related to the clas- 
sical periodic orbits of the system through scmiclassical 
theories like the Gutzwiller's trace formula 6J. Hence, 
it is important to understand the self-affine properties 
of the spectrum since it has implications at the level of 
classical dynamics of the system. 

The idea of characterizing fluctuations and its scal- 
ing is of interest in other areas as well. For instance, 
the surface-height fluctuations in the surface growth 
processes reflect the morphological changes that deter- 
mine the physical and chemical processes such as crystal 
growth, metal deposition etc. [3. Fluctuations in the 
EEG series might tell us about the nature of physiological 
process taking place 8]. In all these cases, fluctuations 
and how it scales provide an important clue to the phys- 
ical process under study. Hence, the results presented 
here have relevance in fields beyond quantum physics. 

In this paper, we take the time series point of view 
for the quantum spectrum and compute the exponent H 
for levels of Lanthanide atoms (Sm, Pm and Nd) and 
RMT ensembles using DFA method 0. Note that the 
value of 7 does not necessarily guarantee self-affine na- 
ture 01 ■ We show that the DFA of order 1 is related 
to A3 statistics of RMT and exploit it to obtain theo- 
retical fluctuation functions and H . We show that the 
time series of level sequences and a random time series 
both with the same value of H have the same fluctuation 
properties. 

We denote a spectrum of discrete levels, belonging to 
either a random matrix or a given atomic system, by 
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Ei,i = 1,2, ... ,n + 1. Then, its integrated level density, 
i.e, the total number of levels below a given level, can 
be decomposed into an average part plus an oscillatory 
contribution, N(E) = N av (E) + N osc (E). Spectral un- 
folding is effected by the transform, Aj = N av (Ei), such 
that the mean level density of the transformed sequence 
{A} is unity. All further analysis is carried out using {A}. 
For instance, the level spacing is, d{ = Aj+i — = 
1, 2, ... ,n. Main object of interest in this paper is the 
fluctuations in the 'time series' of levels given by, 

m 

5 m = ^{di - (d)) = -N osc (E m+ i); m = 1, . . . , n. (1) 

i=i 

Note that S m is formally analogous to a time scries with 
the time index denoted by m |3J. The power spectrum 
S(f) of 8 m displays a power law, i.e, S(f) ~ J -7 , where 

1 < 7 < 2 0- In particular, if the system is regular, 7 = 

2 and for Gaussian ensembles and the chaotic systems, 
7=1. Thus, l/f noise should correspond to H = 0. A 
H = random time series should have a profile similar 
to levels of Gaussian ensembles and atoms as seen in Fig. 
da-c). 

To obtain S m , we needed to unfold the spectrum. This 
puts the spectra from various systems on the same foot- 
ing amenable to direct comparison. For the atomic sys- 
tems, unfolding was done using an empirical function 
to fit the given levels. For the Gaussian ensembles the 
Wigner's semi-circle law 0] was used for unfolding. 

First, we show the results from Lanthanide atoms. The 
atoms in this series exhibit complex configuration mixing 
and spectra. Recently, a series of studies on the statis- 
tical properties of these levels from RMT point of view 
were reported ■ All the calculations for levels of Sm, 
Nd, Pm were done using the GRASP92 code and for de- 
tails we refer to I n Fig- E a ); we show S m for levels 
obtained from Sm. Note that the profile is similar to lev- 
els from Gaussian Orthogonal Ensemble (GOE), shown 
in Fig. ^b), and a random time series with H = ( Fig. 
Hie) ). This is in contrast to the level profile of Gaussian 
Diagonal Ensemble (GDE) ( Fig. Hfd)) defined as an en- 
semble of diagonal matrices whose elements are Gaussian 
distributed random numbers. 

To quantitatively obtain the exponent H, we compute 
it using the detrended fluctuation analysis (DFA) tech- 
nique, which has become a popular tool to study long 
range correlations in time series 9]. DFA is based on 
the idea that if the time series contains nonstationari- 
ties, then the variance of the fluctuations can be studied 
by successively detrending using linear, quadratic, cubic 
... higher order polynomials in a piecewise manner. The 
entire time series is divided into segments of length s and 
in each of them detrending is applied and variance F 2 (s) 
about the trend is computed and averaged over all the 
segments. This exercise is repeated for various values of 
s. For a self-affine fractal, we have, 
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FIG. 1: 5m— m curve for (a) Sm I atomic levels (b) GOE levels 
(c) time series with H = and (d) GDE levels. Note that 
(a,b,c) have qualitatively similar features as opposed GDE 
levels in (d). 



TABLE I: Numerical estimates of Hi and H. 



System 


Hausdorff measure (H) 


m 


Sm 


0.045 


1.045 ± 0.031 


Pm 


-0.025 


0.975 ± 0.024 


Nd 


-0.003 


0.997 ± 0.013 


GOE 


0.0003 


1.0003 ± 0.0083 


GUE 


0.0047 


1.0047 ± 0.0089 


GSE 


0.0025 


1.0025 ± 0.0125 


GDE 


0.5080 


1.5080 ± 0.0044 



F(s) oc s H 



(2) 



where < H < 1. In Tabled we show the numerical 
estimates of H obtained for various atomic and Gaussian 
ensemble cases. The results tabulated alongside repre- 
sent average of F(s) taken over an appropriate ensemble. 
For instance, for the Sm levels, the ensemble consists of 
6 sequences each with 650 levels. Similarly for Pm and 
Nd. In the case of Gaussian ensembles of RMT, the re- 
sults were obtained from ensemble of 50 level sequences 
each with 4600, 3900, 2100 levels respectively for Gaus- 
sian Orthogonal ensemble, Unitary ensemble (GUE) and 
Symplectic ensemble (GSE). For the uncorrelated levels 
from GDE, we get H = 0.508 in good agreement with 
the expected value of 1/2. For the atomic and Gaussian 
ensembles, the value of H is close to zero. In such cases 
standard DFA becomes inaccurate and DFA is performed 
after integrating S m once more so that the new exponent 
is iJi = H +1. Tabled shows that Hi is close to unity for 
atomic and Gaussian ensembles. Thus, the GDE levels 
are self-affine fractal with H = 1/2 but those of atomic 
and Gaussian ensembles (H = 0) do not fall in this class. 

In order to obtain analytical estimate for the value 
of H, we appeal to the ^(s) statistic widely studied 
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in RMT to quantify rigidity of spectrum, i.e, given a 
level, a rigid spectrum will not allow any flexibility in 
the placement of the next level. As(s) is a measure of 
this flexibility or conversely the rigidity. It is defined by, 



A 3 (A, S ) = - 
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i(X) -aX- bfdX (3) 



where n(X) = n av (X) + n osc (X) is the integrated level den- 
sity for the unfolded spectrum {A} and a and b are the 
parameters from linear regression and A is any reference 
level. Now, Aa(s) is an average over various choices of 
A. Notice that n(X) = N(E) and due to unfolding proce- 
dure, n av (X) = A. This implies that, n osc (A) = N osc (E). 
Then, we could rewrite Eq. Q as, 
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A 3 (A, s) = ^- min / [n osc (A) -{a- 1)A - b] 2 dX (4) 

2s a,b J- X _ s 

Now we discretize this equation to relate A 3 (s) with 
the standard fluctuation function F(s) obtained by DFA. 
Substituting from Eq. for n osc (A), we get, 

-, M j+2s 



M ^ 2s 

3=0 m=j + l 
M 

:j-Y^F\s 1 j)=F 2 (s) 

3=0 



(5) 



where Y(X m ) — (a* — l)X m + b* with best fit parameters 
a* and b* and j indexes M different choices of A. Since 
the DFA is applied with linear detrending on segments of 
length s, then Eq. (JSJ) boils down to A 3 (s) = F 2 (s), and 
the subscript denotes DFA(l). Then, ensemble averaging 
and writing -Fi(s) in terms of Aa(s), we have, 



(F 1 ( S )) = (A 3 ( S ) 1 / 2 ) 



(6) 



This main result of the paper shows DFA(l) to be equiva- 
lent to A3 statistic. From RMT, the averages for (A3(s)) 
are well known. For the GDE, we have (Ap^'^(s)) = 
s/15. Then, the result for uncorrelated set of levels is, 



<*i00> 



S \ 1/2 

15/ 



(7) 



On a log-log plot we expect a slope of H = 1/2. This 
is amply verified in Fig. |2td) for GDE with slope 0.508. 
In our numerical simulation, GDE consists of 50 level 
sequences of 15000 levels each. 

The numerical results for atomic levels and Gaussian 
ensembles in Table [Q show that fl w (Hi « 1). To 
make correspondence with H±, we consider modified A3 
statistics given by, 



A(A,s) = - 



1 

mm 

2s a,b 



\+s 



[A n(A) - aX - bfdX 



(8) 



FIG. 2: (Color online) The fluctuation function (Fi 2 (s)} in 
a semi-log plot for (a) atomic levels, (b) Gaussian ensembles 
and (c) random time series with H — and power spectrum 
l/Pf. (d) shows ln(fi 2 (s)} vs Ins for GDE and a random 
walk time series. The slopes of best fit (solid) lines are quoted 
in the text. 



which corresponds to integrating the S m once more before 
performing DFA on it. We obtain the asymptotic random 
matrix average to be [l3[ . 



<A( S )> 



s 2H - cx s 2 



(9) 



for all the three ensembles. Then, using EqEl (Fi (s)) oc s 
and hence the theoretical slope Hi = 1 and corresponds 
to standard Hausdorff measure H — 0. The numerical 
values in last column of Table [I] confirm this RMT based 
result. For Pm and Nd, note that the small negative 
values for H are in the error bar as seen by the error 
estimates provided. 

Using the time series analogy, we explore the nature of 
fluctuations at H = for Gaussian ensembles. This, 
in turn, throws further light on the fluctuation prop- 
erties of H = time series. From RMT we have, 
(A 3 (s)) /3 w i^fF, where j3 = 1,2,4 for GOE, GUE and 
GSE respectively. Substituting this in Eq. ©, the fluc- 
tuation function at H = is, 



(Fi(s)) 



In 2tts 

(3tt 2 



1/2 



(10) 



At H — 0, the ensemble averaged fluctuations display 
logarithmic scaling with length of level sequence s con- 
sidered. For convenience, we plot (Fi 2 (s)} in a semi-log 
plot. Then we expect slopes for level sequences from 
GOE, GUE and GSE to be 1/(3tt 2 . In Fig. 0(a), we 
show results from atomic levels of Sm, Pm and Nd. In 
all the three cases, the straight-line obtained shows that 
the DFA fluctuation function is in good agreement with 
the theoretical result based on As(s) statistics. Their 
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slopes, 0.105, 0.104 and 0.109 for Sm, Pm and Nd agree 
closely with 1/ir 2 . The deviations for Ins > 2.2 occurs 
due to breakdown of universality in the presence of sys- 
tem specific information. This is related to theperiod of 
the short periodic orbits in the atomic system [14J. Thus, 
quantum spectrum of complex atoms can be character- 
ized by H = and display logarithmic scaling. Similar 
results are obtained for the Gaussian ensembles. Fig. 
Elb) shows linear relationship in semi-log plot for levels 
from GOE, GUE and GSE. The measured best fit slopes 
0.0995 (P = 1), 0.0511 (/3 = 2), 0.0257 (J3 = 4) closely 
approximate the theoretical value l/f3ir 2 . The DFA ap- 
plied to spectral levels can distinguish between the three 
ensembles. This is not possible with the power spectrum 
method unless one measures the slope and the intercept 
together. Intercept measurement is often difficult due to 
deviations from scaling for small frequencies. 

What happens to an anti-persistent time series at 
H = ? We generate random time series with H = 
and power spectrum S(f) = l/2j3n 2 f. The technique 
is to generate Gaussian random numbers u and take its 
Fourier transform to obtain u. Then, take the inverse 
Fourier transform of iiy/\S(f)\ to obtain a time series 
with the required properties [15J. By varying (3, we get a 
synthetic random time series that will mimic level fluc- 
tuations of GOE, GUE and GSE. The DFA applied to 
10 member ensemble of such time series (each of length 
5000) is shown in Fig. 0{c). On a semi- log plot, they 
produce the expected linearity and slopes 0.1035, 0.0517 
and 0.0258 agree quite well with the theoretically ex- 
pected result l//?7r 2 . We conclude that an appropriate 
time series with 1// noise is similar to a level spectrum 
from random matrix ensembles and complex atoms and 
both display fluctuations that scale in a logarithmic man- 
ner at H = 0. These results, viewed in conjunction with 
Bohigas et. al. conjecture [lj, mean that we can expect 
the spectrum of classically regular systems to be charac- 
terized by H = 1/2 and chaotic systems by H = 0. 

Even though we motivated our work using quantum 
systems, we emphasize that if the spectrum is viewed as 
a time series, then we can apply our main results to any 



classical system where self-affinity of the time series is 
relevant. For instance, in studies of surface roughness, 
growth models where surface morphology determines the 
properties of the physical and chemical processes. Ex- 
amples of surface growth are the evolution of landscapes, 
vapor deposition, propagation of flame fronts and growth 
of bacterial colonies etc. Q. In these cases, a relevant 
quantity of interest is the correlation function for the 
surface heights h(x), C(r) = (\h(x) - h(x + r)] 2 ) 1 / 2 . In 
many cases, C(r) ~ r H over a large range of length scales, 
i.e, m any such growing surfaces are self-affine fractals. In 
Ref. [l(j a modified form for the correlation function of 
rough surfaces is studied which gives logarithmic fluctu- 
ations, similar to Eq. (|1C)|> obtained above, as H — > 0. 
This has been confirmed by numerical simulations of sur- 
faces as well as from experimental results of thin-film 
evaporation |16| . Another feature of surface growth is 
the existence of roughening transition for certain critical 
parameter values. At roughening transition, the expo- 
nent is H = and C(r) displays logarithmic fluctuations 
as experimentally verified for Ag(115) sample We 
have also shown that H = corresponds to log fluctua- 
tions using random matrix methods. Hence, taking the 
time series perspective and RMT could be beneficial for 
various other problems. 

To summarize, we treat the quantum spectrum like 
a time series. In this picture, the uncorrelated spectral 
levels behave like a random walk series with Hausdorff 
measure H = 1/2 and is a self-affine fractal. For the 
spectrum of complex atoms and Wigner-Dyson random 
matrix ensembles, H = 0, an anti-persistent time series. 
We show that the DFA is equivalent to A3 statistics of 
RMT and we exploit this connection to obtain theoretical 
estimates for H . This shows that at H = the spectral 
fluctuations display logarithmic behavior, a feature seen 
in many experimental and model simulations in surface 
growth studies. Thus, H = is an unusual and interest- 
ing limit and will require further investigation. 

We thank Prof. J. C. Parikh, Dr. P. Panigrahi, Prof. 
V. B. Shcorey and P. Manimaran for useful discussions. 
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